Experimental Violation of Two-Party Leggett-Garg Inequalities 
with Semi-weak Measurements 
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We generalize the derivation of Leggett-Garg inequalities to systematically treat a larger class of 
experimental situations by allowing multi-particle correlations, invasive detection, and ambiguous 
detector results. Furthermore, we show how many such inequalities may be tested simultaneously 
with a single setup. As a proof of principle, we violate several such two-particle inequalities with data 
obtained from a polarization-entangled biphoton state and a semi-weak polarization measurement 
based on Fresnel reflection. We also point out a non-trivial connection between specific two-party 
Leggett-Garg inequality violations and convex sums of strange weak values. 

PACS numbers: 42.50.Xa,03.65.Ta,42.50.Dv 



To better understand and identify the apparent di- 
vision between macroscopic and microscopic behavior, 
Leggett and Garg have distilled common implicit as- 
sumptions about the macroscopic world into a set of ex- 
plicit postulates that they dub macrorealism (MR) 
From these postulates, they construct inequalities anal- 
ogous to Bell inequalities [2| but involving multiple cor- 
relations in time. Such Leggett-Garg inequalities (LGls) 
must be satisfied by any theory compatible with MR, but 
may be violated by quantum mechanics. As such, LGl 
violations have received increasing interest as signatures 
of distinctly quantum behavior in qubit implementations 
[^-Q, and have been recently confirmed experimentally 
in both solid-state iG] and optical systems [7[. 

In this Letter, we demonstrate a technique for system- 
atically deriving generalized LGls that admit multiple 
parties, invasive detection, and/or ambiguous detector 
results by considering a specific two-particle experimen- 
tal setup with three measurements. We proceed to exper- 
imentally violate several such two-party LGls simultane- 
ously with a single data set produced from a setup using 
a semi-weak polarization measurement on an entangled 
biphoton state. The contextual values (GV) analysis of 
quantum measurement Q suggests a direct comparison 
between the classical and quantum treatments. Finally, 
we show that specific two-party LGls are equivalent to 
constraints on convex sums of conditioned averages (GA), 
which are the generalizations of the quantum weak value 
to an arbitrary measurement setup d, Q . The technique 
may be easily extended to check data from a setup with 
any number of measurements and parties. 

Generalized LGls. — A MR theory consists of three key 
postulates: (i) if an object has several distinguishable 
states available to it, then at any given time it is in only 
one of those states; (ii) one can in principle determine 
which state it is in without disturbing that state or its 
subsequent dynamics; and, (iii) its future state is deter- 
mined causally by prior events [ij. Furthermore, we ac- 
knowledge that physical detectors may be imperfect by 
being (a) invasive by altering the object state during the 
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FIG. 1. MR measurement schematic. An object pair is picked 
from an ensemble ( at time to. At ti object 1 of the pair in- 
teracts with an imperfect detector for the property Ai, which 
reports a generalized value ai. At t2 both objects interact 
with unambiguous detectors for the properties Bi and B2 that 
report values bi and 62. The two-party LG correlation C is 
constructed from the measured results. 

interaction, or (b) ambiguous by reporting results that 
only correlate probabilistically with the object state due 
to inherent detector inefficiencies or errors. 

For convenience we consider dichotomic properties in 
what follows, though the discussion can be easily ex- 
tended. Unambiguous detector outcomes will be assigned 
the (arbitrary) values {—1,1} corresponding to the two 
possible states of the property being measured. Ambigu- 
ous detectors will be calibrated to report the same en- 
semble average as an unambiguous detector for the same 
property. To do so, their outcomes must be assigned 
generalized values a € S from an expanded set S, with 
min S < —1 and max 5* > 1 , to compensate for the im- 
perfect state correlation of the outcomes. Such general- 
ized values are the classical equivalent of quantum CV Q 
and may be determined by measuring pure ensembles of 
either ±1. 

We now derive a specific two-party generalized LGl for 
a particular experimental setup, keeping in mind that the 
method may be extended to any setup. Consider a pair 
of MR objects that interacts with a sequence of detec- 
tors as shown in Fig. [1] At time to the pair is picked 
from a known ensemble ^. At time ti object 1 of the pair 
interacts with an imperfect detector for the dichotomic 
property Ai, which reports a generalized value ai € Si. 
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Finally, at time t2 objects 1 and 2 interact with unam- 
biguous detectors for the dichotomic properties Bi and 
B2, respectively, which report the values &i, &2 <= {^l) !}■ 

For each object pair, we can keep all three results 
to construct the correlation product 016162, or we can 
ignore some results as non-selective measurements 
to construct the alternate quantities ai, 61, 62, ai6i, 
q;i62, or 6162. Since the latter terms involve voluntary 
loss of information after the measurement has been per- 
formed, we can compute them all from the same data 
set. Exploiting this freedom, we construct the correlation 
C = ai + aibib2 — 6162 for each measured pair, which 
lies in the range, — 11 — 2minS'i| < C < 1 2 max 6*1 — 1|. 

We repeat this procedure many times and 
average the results of C to obtain, (C) = 
Ec.MM Piai\0Pibub2\C, ai) (ai + aib,b2 - 6162), 
where P(ai|C) is the probability of detecting ai given 
the initial ensemble C, and P(6i,62|C,ai) is the proba- 
bility of detecting 61 and 62 given the initial ensemble C 
and the possibly invasive detection of ai. 

Generally, we cannot separate the sums due to the 
Q!i-dcpcndcncc of P(6i, 62IC, ai), so the best guaranteed 
bounds are — 1 1 — 2 min Si\ < (C) < 1 2 max 5*1 — 1 1 . As a 
special case, if the detector for ^1 is unambiguous then 
min5i = — 1, maxS'i = 1, and wc find the LGI, 

-3<{Ai+AiBiB2-BiB2)<l. (1) 

Alternatively, if the detector is noninvasive so that 
-P(6i, 62IC, ai) — P(6i,62|C) then the sums do sep- 
arate and we can average Ai first to find, (C) = 
Ebi,6,m,&2|C)((^i)(l + ^i&2)-6i62). Since -1 < 
(^1) < 1, each term can take only three possible val- 
ues {—3,-1, 1} and we again recover ([T]). Therefore, any 
violation of ([T|) will imply that at least one of the postu- 
lates (i-iii) of MR does not hold, or that the detector for 
Ai is both invasive and ambiguous. 

We can construct many such LGIs from the same data. 
For example, the three detectors in Fig. [T] allow the con- 
struction of the 2"^ — 1 nontrivial correlation terms listed 
earlier, which can be combined with the three coeffi- 
cients { — 1,0,1} Ignoring an overall sign, wc can 
construct (3^ ~^ — l)/2 = 1093 nonzero LGI correla- 
tions bounded in a similar manner to ([1]). The subset 
of (32 -1 - l)/2 = 13 single-object LGIs can be obtained 
by ignoring the B2 detector. Furthermore, if a fourth 
detector for A2 were added before the detector for B2, 
we could test (3^*-^ - l)/2 = 7174453 such LGIs. One 
is formally identical to the CHSH-Bell inequality (see 
also [13 ): but tests MR and not Bell- locality. 

For contrast, the original approach in [l[ combines sep- 
arate experiments for each correlation between ideal de- 
tectors to form a single LGI. Our approach uses a single 
experimental setup to determine all 2*^ — 1 correlations 
between M general detectors to form a large number of 
LGIs. Hence we obtain an exponential improvement in 
experimental complexity for large M. 
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FIG. 2. (color online) Experimental setup. A 488 nm laser 
produces degenerate down-converted photon pairs. The po- 
larization of the photon in the lower arm is rotated by 45° 
with a half-wave plate, then undergoes semi-weak polarization 
measurement in the {h,v} basis using Fresnel reflection (Ai) 
that encodes the information in the resulting spatial modes, 
and is finally projected into the {6,9±} basis with polariz- 
ers set at angle 6 (Bi). The polarization of the photon in 
the upper arm is projected into the {h, v} basis with another 
polarizer (i?2). The half and quarter waveplates prior to the 
polarizers are used for tomography of the input state; during 
data collection they are removed from the lower arm and used 
to switch between h and v polarization in the upper arm. 

Conditioned averages. — A single-object LGI, —3 < 
{Ai + AiBi — Bi) < 1, was considered in Q and shown to 
have a one-to-one correspondence with an upper bound 
to the average of Ai conditioned on the positive value of 
Bi'. i{A) < 1. Three other LGIs similarly correspond to 
the bounds i{A) > — 1, and — 1 < -i(^) < 1, as checked 
experimentally in~|3]. 

Wc now extend these results to the two-object case 
using ID). First we define a marginal probability 
of measuring 61 and 62 given any result of Ai as 
P(6i,62|C,^i) - Ea,Piai\0Pibub2\C,ai). Then we 
define a conditional probability of measuring ai given 
the measurement of 61 and 62 as, P(ai |C, 61, 62) = 
P(ai|C)P(6i,62|C,ai)/-P(6i,62|C,^i). Therefore, the av- 
erage of Ai conditioned on the measurements of 61 and 
62 is 61,62(^1) = Eai P{ai\C,bi,b2)ai. 

Using this definition, we rewrite the upper bound of ([1]) 

as Ehi,6,^(fel,^2|C,^l)(fcl,62(^l)(l+fel&2)-6i62) < 1 

and insert the possible values for 61 and 62 to find the 
CA constraint, 

i,i(^i)p++ _i,_i(^i)p- < 1, (2) 

where p± = P(±l, ±1)/(P(1, 1) -I- P(-l,-l)), and 
P{hi) — PihiKi-^i)- The degeneracy of the product 
value 6162 results in an upper bound for a convex sum of 
CAs, in contrast to the single-object result in j^]. A suffi- 
cient condition for violating ^ is for both CAs to exceed 
1 simultaneously. Conversely, if all CAs were bounded by 
1, then it would be impossible to violate ([2]) or ((T|). 

Quantum formulation. — Projective quantum measure- 
ments produce averages of eigenvalues analogous to the 
results of an unambiguous detector, but non-projective 
quantum measurements produce averages of contextual 
values 3 which need not lie in the eigenvalue range and 
arc therefore analogous to the results of an ambiguous 
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FIG. 3. (color online) Real (left) and imaginary (right) parts 
of the reconstructed density matrix in the {h, v} basis. Yellow 
and red represent positive and negative values, respectively. 

detector. By measuring Ai weakly we can find quantum 
mechanical violations of ([T]) and 

Specifically, if we start with a 2-object density op- 
erator p and measure Ai generally such that Ai = 
Sai «inai = "i-E'qi (where {ai} are the eigenval- 
ues corresponding to the projections {Hai} and {ai\ are 
the CV corresponding to the POVM = Ml^M^i}), 
and then measure B1B2 projectively such that Bi®B2 = 
62 ^i^sll&i ® flfc^, we will find that the average cor- 
relation (C) = (^1 + A1B1B2 - B1B2) has the form, 

(C)= Piai;bub2\p)iai+aibib2-hb2), (3) 



ai ,bi ,62 



where P(q;i; &i, 62I/0) = Tr 



P 



the probability of measuring outcome ai of the general 
measurement of A^ followed by a joint projection of &i&2- 
The appearance of the CV instead of the eigenvalues of 
A in combined with the non-separable probability 
P{ai; hi, &2I/5) allows violations of the LGI ([IJ. 

The left side of ([2]) follows from ([3]), where 
P{bi,b2\p,Ai) = Y.o.,P{(^i-MM\p) and 61,62(^1) = 
X;^! ai^'(ai;&ij^2|/5)/^'(^^i,&2|p, ^1) is a quantum CA 
as defined in that converges to a weak value 0] in 
the limit of minimal measurement disturbance. 

Experimental setup. — To implement Fig. [1] we use the 
polarization of an entangled biphoton with the setup 
shown in Fig. [21 A glass microscope coverslip mea- 
sures a Stokes observable Ai semi- weakly as described 
below, and polarizers measure Stokes observables Bi and 
B2 projectively. We produce degenerate non-colinear 
type-II down-conversion by pumping a 2 mm walkoff- 
compensated BBO crystal [13| with a narrowband 488 
nm laser. The down-converted light passes through auto- 
mated polarization analyzers and 3 nm bandpass filters at 
976 nm in each arm before being coupled into multimode 
fibers connected to single photon avalanche photodiodes 
(SPAD). We detect coincidences using a 3 ns window. 
We perform state tomography with maximum likelihood 
estimation which gies the state shown in Fig. |3] with 
concurrence C = 0.794, and purity Tr [p^] = 0.815, and 
which resembles the pure state j?/;) = (1^"^) +i\vh))/^/2. 

After the state tomography, we remove the half- and 
quarter-wave plates from the lower arm and insert either 
a mirror or a coverslip using a computer-controlled trans- 




FIG. 4. (color online) In all data plots, solid lines indicate the- 
ory and points indicate experimental data, a) (iji^') (green, 
fiat), e((Ti^') (red, decreasing), and (o-i""^') (blue, increas- 
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b) e,h{(Ti^') (red, bottom right), and e^,t,(i5-i"'^') (blue. 



bottom left), violating negative bounds, unlike 9^,h{(Ti^') (or- 
ange, increasing), and (green, decreasing). 

lation stage. The reflected light passes though a polar- 
ization analyzer and couples into a third fiber and SPAD. 
We align the coverslip and the mirror to be parallel with 
an incidence angle of 40° relative to the incoming beam. 
Finally, we optimize the fiber incoupling and balance the 
collection efficiencies with attenuators so that the coin- 
cidences between the upper arm and either of the lower 
arms differ by only a few percent when the mirror is taken 
in and out of the beam path. 

The coverslip acts as a polarization-dependent beam- 
splitter measuring Ai = cfz- Averaging over the 3 nm 
bandwidth and the thickness variation 150 ± 0.6 /im) 
produces an average Fresnel reflection similar to that of 
a single interface, with horizontal ih) polarization rela- 
tive to the table exhibiting zero reflection near Brewster's 
angle and vertical (i>) polarization exhibiting increasing 
reflection with incident angle. 

For a pure state of polarization ~ a\h) + I3\v) with 
jap -I- = 1, the resulting state after passing through 
the coverslip is \ip') = {loi\h) + fil3\v))\r) — {'^oi\h) + 
irjj3\v))\t), where j € {tif}i specify the transmit- 
ted and reflected spatial modes of the coverslip, and the 
reflection and transmission probabilities for h- and v- 
polarized light are Rh = 7^, Rv = Tf , = 7^, and 
Tv = ?7^, such that Ri + Ti = 1. Written this way, the 
coverslip reflection can be viewed as a generalization of 
the weak measurement in [l5| and discussed in Q . 

From \ip'), we find the measurement operators for 
the back-action of the coverslip outcomes to be Mr = 
^Ilh+fjlly and Mt ~ 711/1-1-77111,, where Hi, i G {/i, u}, are 
polarization projectors. The corresponding POVM ele- 
ments are Er = Rh^h + Rv^v and Et = Tfjlh + 71>n„, 
with which we can expand the polarization Stokes op- 
erator as (T2 = Hh — Hv = ttr-Er + atEt, as discussed 
before ([3]), where a,. = {T^ + T^)l{Rh — Rv) and at = 
-{Rh + Rv)/{Rh - Rv) are the CV. 

We determine the values of Rh and Rv with cali- 
bration polarizers before the coverslip, yielding Rh = 
0.0390 ± 0.0007 and Rv = 0.175 ± 0.001. The reflected 
arm is largely projected to w, while the transmitted arm is 
only weakly perturbed, making the total coverslip effect 
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FIG. 5. (color online) LGI correlation {—ai^^ — o"i^V^^''(Ti^' — 
Cg^Vi'^') (red, squares) and the corresponding convex sum of 

the CAs 9,h{—S'i^'') and ej_,v{—o'i^^) (blue, circles), both vio- 
lating their upper bounds of 1 in the same domain of 6. Com- 
pare to Fig. |4] (b) and note that the LGI violation includes 
the region where the two CAs both exceed their bounds. 

a semi-weak measurement. The CV, = —13.1 ± 0.1 
and at = 1.57 ±0.01, are correspondingly amplified from 
the eigenvalues of CTz. 

Results. — To complete the state preparation, we place 
a half-wave plate before the coverslip in the lower arm 
and rotate the polarization by 45° to produce a state 
similar to lip") = {\ha) + i\vd})/V2. We then measure 
([1]) by choosing the observables Ai, Bi, and B2 to be 
the Stokes observables and (tI^^ , respectively, 

where as is the CTz operator rotated to the {9,9±} basis 
(e.g. ao° = o'z and 1T450 = ax). By changing the single 
parameter, 9, we can explore a range of observables. 

Fig. |4]shows the various averages of (ji^-* . Averaging all 
results for orthogonal settings on a^^^ and ai^^ gives the 

expectation value ((ji^^), which is properly constant and 
near zero for all 9 since the reduced density operator is 
almost fully mixed. Averaging only the results for the or- 
thogonal settings of a'i^'' gives the sin glc CAs g{ai^^) and 
which are also well-behaved. Finally, averaging 
only the results for specific settings gives the double CAs 
9,v{^^/^), 9_L,h{^P), 0,-u('5'i"^^), and e^.tiio-i^^), which can 
exceed the eigenvalue range for some range of 9 due to the 
non-local correlations in the entangled biphoton state. 

Using the same set of data. Fig. [S] shows the upper 
boundoftheLGI-3< (-ai'^-txi'V^'Vi'^-^'Vi'') < 
1 being violated in the same range of 9 that the appro- 
priate convex sum of CTz^^) and dj^^vi^ai^"^) violates 
its upper bound according to ([2|) . 

We can violate several more LGIs using the same set 
of data as well. Fig. [S] shows two such correlations, 
(ai^Vi^)+ai^V«-ai^VW), and {-ai'^ ai'^ a^'^ + 
that between them violate an upper bound over 
nearly the whole range of 9, for illustration. 

All solid curves in Figures |4l[5l and [6] are quantum pre- 
dictions analogous to ([3]) using the measurement opera- 
tors, CV, and the reconstructed initial state. They also 
include compensation for a few percent deviation in the 
thickness of the half-wave plate in the upper arm. The 



FIG. 6. (color online) LGI correlations {ai^^ai^^ +o-i^^ag^^ — 
cri'V*'') (red, circles), and (-ai'Vi'^ + ai^V^'' + ffi'V*'') 
(blue, squares) violating their upper bounds of 1 for nearly 
the entire 9 domain. 

points indicate experimental data and include Poissonian 
error bars. The small discrepancies between theory and 
data can be explained by sensitivity to the state recon- 
struction and additional equipment imperfections. The 
violations indicate either that MR is inconsistent with 
experiment or that the semi- weak measurement device is 
both invasive and ambiguous in the MR sense. 

Conclusion. — We have illustrated the derivation of 
generalized single-setup LGIs allowing for multiple par- 
ticles and measurements with more realistic detectors 
by considering a two-particle example, and have demon- 
strated simultaneous violations of several such two-party 
LGIs using the same data set from a biphoton polariza- 
tion experiment. Due to the single setup, any dataset 
may be similarly examined for inherent LGI violations. 
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